Molecular stability of a vaccine strain of Canine coronavirus 
after serial passages in A72 cells 


Estabilidade molecular de uma amostras vacinal de Coronavirus 
canino apos passagens seriadas em células A72 


Iracema Nunes de BARROS!”; Sheila Oliveira de Souza SILVA!”; Sueli Akemi TANIWAKI?”; 
Paulo Eduardo BRANDAO!” 


"Universidade de Sao Paulo, Faculdade de Medicina Veterinaria e Zootecnia, Departamento de Medicina Veterinaria Preventiva 
e Satide Animal, S40 Paulo — SP, Brazil 
? Universidade de Sao Paulo, Coronavirus Research Group, Sao Paulo - SP, Brazil 


Abstract 

Canine coronavirus (CCoV) exists in types I and II and infects dogs leading mainly to enteritis, though type II has 
already been associated with generalized and highly lethal infection. A CCoV-type II inactivated vaccine produced in 
A72 canine cells is available worldwide and largely used, though the molecular stability after serial passages of vaccine 
seeds is unknown. This article reports the evolution of the CCoV-II vaccine strain 1-71 in A72 cells based on partial 
S gene sequencing, showing the predominance of neutral evolution and the occurrence of four sites under purifying 
selection. Thus, cell-adapted strains of CCoV-II may be genetically stable after serial passages in a same cell line due to 
a stable virus-host relationship. 
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Resumo 

O Coronavirus canino (CCoV) ocorre como tipos I e II e infecta caes, levando principalmente a enterite, apesar do 
tipo II ja ter sido associado a infeccao generalizada e altamente letal. Uma vacina de CCoV-II inativada produzida 
em células caninas A72 é disponivel mundialmente e largamente utilizada, apesar da sua estabilidade molecular apés 
passagens seriadas de sementes vacinais ser desconhecida. Este artigo relata a evolucgao da amostra vacinal CCoC-II 
1-71 em células A-72 com base em sequenciamento parcial do gene S$, demonstrando predominio de evolugao neutra 
e a ocorréncia de quaro sitios sob selegao purificante. Portanto, amostras de CCoV-II adaptadas a cultivos celulares 
podem ser estaveis geneticamente apos passagens seriadas em uma mesma linhagem celular devido a existéncia de uma 
relacdo estavel virus-hospedeiro. 
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The S protein is the main target of neutralizing antibodies 
and the most polymorphic protein in coronaviruses. It 
occurs as trimmers, with a highly variable ectodomain 
(subunit $1) and a more stable endodomain (S2 subunit) 
involved in membrane fusion (BOSCH et al., 2003; 
MASTERS, 2006). 

Estimates on the time of the most recent common 
ancestor (tMRCA) place the origin of Alphacoronaviruses 
293 million of years ago and, due to the proofreading 
activity of the NSP14 of coronaviruses, these viruses show 
a number of mutations per site per replication up to E-5, 
which is quite similar to the one found in DNA viruses 
(WERTHEIM et al., 2013). 

In spite of the large use of cell-passaged CCoV vaccine 
seeds for years on and the use of high-passage of vaccine 
strains in vaccine trials (PRATELLI, 2007), no data on 
the molecular diversity and the mechanisms of molecular 
evolution of CCoV in vitro is available. This knowledge is 
paramount for the understanding of the genetic stability 
of vaccine seeds and might also help shed some light on 
basic mechanisms of coronavirus molecular evolution. 

The aim of this investigation was to assess the clonal 
diversity of a CCoV-II vaccine strain after serial passages 
in A-72 cells based on the molecular diversity of the $2 
subunit of the spike gene and the selection regime driving 
CCoV in vitro evolution. 


Material and Methods 


CCoV strain and cell line 

CCoV 1-71 type II vaccine strain (ATCC VR 809, 
10*°TCID,,/mL) and canine fibroma A72 (ATCC CRL 
1542) cell line were used for this study and were kindly 
provided by Biovet Laboratories, Brazil. 

Cells were propagated in Eagles MEM supplemented 
with 10% fetal bovine serum (FBS) and 25mM HEPES 
(4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid). 


Serial passages in A72 

For the first passage, ImL of CCV 1-71 strain diluted 
1:10 (v/v) in Eagle’s MEM was inoculated in 25 cm? flasks 
containing partially (~70%) confluent A72 cells and 
incubated at 37°C/ 1 hour for virus attachment. Next, the 
virus suspension was discarded and 7ml of Eagle’s MEM 
supplemented with 2% FBS was added to the flask. 

Cells were daily observed with a light microscope for 
cytopathic effect (CPE) for five days. After this observation 
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period, the monolayers were frozen at -80°C regardless 
of the presence of CPE. Then, 1mL was used in the next 
passage following the same protocol described above up 
to the 9" passage 

For each passage, a 25 cm? with A72 cell was mock- 
infected (negative control) with Eagle's MEM without virus 


under the same conditions described for virus inoculation. 


RT-PCR for partial S gene amplification 

For the original CCoV 1-71 and each passage (CCoV 
and negative control), total RNA was extracted with TRIzol 
Reagent™ (Life Technologies) following the manufacturer’s 
instructions. cDNA synthesis was carried out with Random 
Primers” (Life Technologies ) and M-MLV Reverse 
Transcriptase” (Life Technologies) accordingly to the 
manufacturer’s instructions. 

For partial S gene amplification, the primers 
described by Pratelli et al. (PRATELLI et al. 
2004) (S 5S TGCATTTGTGTCTCAGACTT37/  S6 
5SCCAAGGCCATTTTACATAAG3’ nt 3492 to 4185 of 
S gene, S2 subunit in the reference strain CCoV 1-71 
GenBank JQ404409) were used with Platinum Taq 
Polymerase High Fidelity” (Life Technologies) following 


the manufacturer’s instructions. 


DNA cloning and Sequencing 

The 694 bp amplicons were purified using ExoSAP- 
IT™ (Affymetrix) and the total DNA quantified with Low 
Mass DNA Ladder™ (Life Technologies) according to the 
manufacturer's instructions. 

The purified amplicons were inserted into the plasmid 
vector pTZ57R/T™ (InsTAclone) in an insert: vector proportion 
of 3:1 and the ligation reaction was performed at 15°C for 16 
hours, followed by the addition of 5 uL of ligation reaction to 
100 uL of competent DH5a E. coli cells grown in LB (Luria- 
Bertani) liquid medium and incubated on ice from 20 to 30 
min. After a thermal shock (42°C/2 min and ice/2 min), 900 uL 
of LB medium without antibiotics were added to the reaction 
and incubated at 37°C / 150 rpm / 1h30 min. 

The cells were then plated on solid LB medium 
containing 100 ug/ml ampicillin, IPTG (5-bromo-4-chloro- 
3-indolyl-B-D-galactosidase), and X-Gal (isopropylthio-B- 
Dgalactoside), and then incubated for 16 hours at 37°C. 

Individual white colonies, predicted to have 
inserts, were selected and subjected to a PCR reaction 
to confirm the insertion using GoTaq Green Master 


Mix” (Promega) accordingly to the manufacturer’s 
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instructions and the  plasmid-targeting primers 
M13F 5S'GTTTTCCCAGTCACGAC3*/ M13 
5‘CAGGAAACAGCTATGAC 3’). 

A number from 9 to 13 clones for each passage and 
of the original virus confirmed as carrying the inserted 
amplicon were subjected to DNA sequencing with BigDye 
3.1™ (Applied Byosystems) and ABI-3500 Genetic Analyzer™ 
(Applied Byosystems) as per manufacturer's instructions, 
using primers M13F and M13R. Evaluation of the quality of 
chromatograms obtained for each DNA strand sequences was 
performed using Phred application online (http://asparagin. 
cenargen.embrapa.br/phph/) and manual checked using the 
program Finch TV (°Geospiza) to search for interpretation 
errors and discrepancies between each DNA strand 
sequenced. The final sequence of each clone was obtained 
using Cap-Contig, with PHRED score > 20, in the program 
Bioedit version 7.2.5 and submitted to BLASTn in http:// 
www.ncbi.nlm.nih.gov/BLAST for homology checking. 


Mutation analysis 

The occurrence of recombination was assessed with 
RDP, GENECONV, Bootscan, MaxiChi, Chimaera, SiScan, 
and 3Seq methods with Bonferroni correction and highest 
acceptable P-value of 0.05 using RDP 4 836 (MARTIN et 
al., 2010). 

Sequences of clones of each passage and the original 
virus (e.g. intra-passage analysis) were aligned using 
Clustal/W in Bioedit version 7.2.5 and compared in terms 
of nucleotide and amino acid identities. Identical sequences 
that predominated in the intra-passage comparison were 
named as Master Sequence for that passage, while those 
in lower frequencies were named Variants. 

Next, the Variants and Master Sequences for all passages 
and the original virus were aligned together including 
Alphacoronavirus homologous sequences retrieved from 
GenBank (JQ404409; JF682842; GQ477367; KC175340; 
JQ404410; EU924790; AY342160; D13096; EU924791; 
EU856361; EU856362; DQ201447; HQ462571; GQ152141; 
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NC_002306; DQ811789; DQ112226) using CLUSTAL/W 
running in Bioedit 7.2.5. 

All sequences generated in this study were deposited 
in the GenBank under the accession numbers: KP281485 
to KP281596. 


Selection regime analysis 

For each codon, estimates of the numbers of 
inferred synonymous (s) and nonsynonymous (n) 
substitutions and the numbers of sites that are estimated 
to be synonymous (S) and nonsyonymous (N) were 
produced using the joint Maximum Likelihood (ML) 
reconstructions of ancestral states under a Muse-Gaut 
model (MUSE; GAUT, 1994) of codon substitution 
and General Time Reversible model of nucleotide 
substitution. For estimating ML values, a tree topology 
was automatically computed. The test statistic dN - dS 
was used for detecting codons that have undergone 
positive selection, where dS is the number of 
synonymous substitutions per site and dN is the number 
of nonsynonymous substitutions per site. Values of p 
less than 0.05 were considered significant at a 5% level. 
These analyses were conducted using HyPhy software 
v. 2.2.1 (POND et al., 2005) implemented in MEGA6 
(TAMURA et al., 2013). 


Results 


Mutation analysis 

A total of 112 viable DNA sequences were obtained 
from passages 1 to 9 and from the original CCoV 1-71 
strain, with a number of clones per passage ranging from 
9 to 13 and only 12 nucleotide mutations found amongst 
these sequences, ranging from 0 to 3 per passage (Table 1) 
but in a passage-exclusive way. No recombination was 
found amongst these sequences. 

From these mutations, five were silent transitions 
(A3519G, T3813C, C3699T, T4155C, and A3527G); 
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three were missense transitions (A3817G, G3868A, and 
G3875A), one silent transversion (A4165C) and two 
missense transversions (G3528T and T3803G). 

Besides, a deletion (TTGATATTA/ Ile/Asp/lle, 
nucleotides 3935 to 3943) was found in one out of 11 clones 
of passage 3. 
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Selection regime analysis 

None of the 12 sites with mutations (Table 1) was 
under positive selection and only sites A3519G, T3813C, 
C3699T, and T4155C were found to be under purifying 
selection, while the remaining eight mutations were 


neutral, including the deletion. 


Table 1 - Nucleotide and putative amino acids mutations of the spike gene of CCoV stain 1-71 passaged in A72 cells. P- passage 
number (0 refers to the original virus); NC - number of clones sequenced per each passage; M/P - number of mutations 
per passage regarding the Master Sequence. Numbers in cells of each mutation represent the number of occurrences of 
that specific mutation. Positions at S gene in the reference strain 1-71 GenBank JQ404409 — Sao Paulo - 2015 


p NC mp A3519G A3527G G3528T C3699T 1T3803G 1T3813C A3817G G3868A G3875A T4155C A4165C sub wasn 
Arg/Arg Glu/Gly Glu/Asp Thr/Thr Phe/Cys Thr/Thr Arg/Gly Cys/Tyr Glu/Lys Asn/Asn_ Thr/Thr Deletion 

(Ile+Asp+lle) 

0 9 2 1 0 0 0 0 1 0 0 0 0 0 0 

1 10 3 0 0 1 1 0 0 0 0 0 0 1 0 

2 11 0 0 0 0 0 0 0 0 0 0 0 0 0 

3 11 2 0 0 0 0 0 0 0 0 0 1 0 1 

4 10 1 0 0 0 0 1 0 0 0 0 0 0 0 

5 13 1 0 0 0 0 0 0 1 0 0 0 0 0 

6 13 2 0 0 0 0 0 0 0 1 1 0 0 0 

7 12 0 0 0 0 0 0 0 0 0 0 0 0 0 

8 11 0 0 0 0 0 0 0 0 0 0 0 0 0 

9 12 1 0 1 0 0 0 0 0 0 0 0 0 0 
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Figure 1 - Clones’ sequences with amino acids mutations in Alphacoronavirus-1 $2 subunit (aa 1172 to 1388 in strain 1-71 Genbank 
JQ404409) after serial passages of strain 1-71 in A72 cells 


Source: this article 
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Discussion 

After nine serial passages in A72 cells, the 1-71 vaccine 
strain of CCoV showed a low relative number of mutations 
when the clonal diversity of a segment of the $2 subunit 
was focused, without a passage-association pattern, which 
might indicate that this strain have reached a quite stable 
virus-cell relationship with no further mutations being 
fixed even after serial passages. 

Though the region of the S gene targeted in this 
study is responsible for membrane fusion and syncytial 
formation and has a higher mutation pressure than the 
S1 region, responsible for binding to cellular receptors 
(LAI; CAVANAGH, 1997), nucleotide polymorphism may 
also be found in $2 region (GALLAGHER; BUCHMEIER, 
2001) allowing also virus entry into a variety of types of 
cells in trypsin-independent way (BELOUZARD et al., 
2009; BORUCKT et al., 2013). 

A study using a FIPV (Feline infectious peritonitis 
virus, a cat alphacoronavirus) strain passaged 100 times in 
A72 cells has shown that, from the original strain to the 8" 
passage, amino acids substitutions occurred in $1 (T599A 
and A743V), while in $2 an amino acid substitution occurred 
only at the 100" passage (R1325E) (PHILLIPS et al., 2013). 

TGEV (Transmissible gastroenteritis virus), another 
alphacoronavirus, has been shown to have a $585A in 
S1 at passage 50 in PK (Porcine Kidney) cells leading to 
attenuation of virulence (LI et al., 2010). 

As seen in Table 1, the mutations detected amongst 
the passages were not fixed, which could be due either to 
genetic drift or to a bias on the randomness of the selection 
of clones for DNA sequencing. 

Though _ these 
alphacoronaviruses are genetically quite stable after a high 


reports provide evidence that 
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CCoV might have been a result of this stability mainly 
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In summary, cell-adapted CCoV-II vaccine strains 
may be highly genetically stable after serial passages 
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